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Abstract. We consider Apollonian circle packings of a half Euclidean plane. We give necessary 
and sufficient conditions for two such packings to be related by a Euclidean similarity (that is, by 
translations, reflections, rotations and dilations) and describe explicitly the group of self-similarities 
of a given packing. We observe that packings with a non-trivial self-similarity correspond to positive 
real numbers that are the roots of quadratic polynomials with rational coefficients. This is reflected 
in a close connection between Apollonian circle packings and continued fractions which allows us 
to completely classify such packings up to similarity. 
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1. Introduction 

A circle packing in M 2 is a set of circles in the plane whose interiors (suitably interpreted) are 
mutually disjoint. An Apollonian circle packing V has the property that for any three mutually 
tangent circles in V, the two circles in the plane that are tangent to all three of them also lie in V . 
Note that our notion of circle includes straight lines where we consider parallel lines to be tangent 
at infinity. 

These types of circle packings have been extensively studied by Graham, et al [21 IH [6] , with a 
focus on those packings for which all the circles have integer curvatures. 

There are four basic shapes that an Apollonian packing may take, and these are illustrated in 
Figure [T] A bounded Apollonian packing (Figure 1(a) ) is a packing V for which a single circle in 



V bounds the entire packing. Here the 'interior' of the bounding circle is the unbounded component 
of its complement. 



A half-plane packing (Figure 1(b) ) is an Apollonian packing V for which at least one of its circles 
is a straight line. The line partitions the plane into two half-planes: one is 'packed' by V, while 
the other is the 'interior' of the line. A special type of half-plane packing is a strip packing 



(Figure 1(c)), in which two of the circles are (necessarily parallel) lines and the remaining circles 



lie in the strip between them. 



An unbounded packing (Figure 1(d) ) is an Apollonian packing which contains no bounding circle 
and no straight line. 

Stereographic projection allows us to relate circle packings in the plane to those on the sphere. 
The four possible configurations in Figure [T] correspond to projection from (a) an interior point, 
(b) a point on only one circle, (c) a tangency point, and (d) a point not on or inside any circle, 
respectively. 

In this paper, we consider the similarity relation on Apollonian circle packings. A similarity is 
a transformation of the Euclidean plane that preserves ratios of lengths. Such a transformation is 
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(c) Strip (d) Unbounded 

Figure 1. Different boundedness properties for Apollonian packings. 



necessarily a composite of a translation, rotation, reflection and/or dilation. Two packings V and 
V' are similar if there is a similarity of the plane that takes circles in V bijectively to circles in V' . 
There may be non-trivial similarities from V to itself, in which case we say that V is self-similar 
and we consider its group of self-similarities. 

Our results concern only half-plane packings. We give a necessary and sufficient condition for two 
such packings to be similar, a classification of the self-similar packings, and a description of all of 
the self-similarity groups. 

We also answer the more specific question of whether two packings are similar via an orientation- 
preserving similarity (that is, one with positive determinant) or via an orientation-reversing 
similarity. Our classification tells us which packings possess an orientation-reversing self-similarity. 

We observe that any half-plane packing V is similar to a packing containing three circles in the 
configuration shown in Figure [2J where L is the x-axis, and a 2 and 1 refer to the curvatures of the 
circles they label. For a > 0, we define V a to be the unique Apollonian circle packing containing 
that configuration. Because each half-plane packing V is similar to such a packing, we restrict 
our attention to studying the packings V a , and we state results in terms of this particular class of 
packings. 
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Our first result relates similarities of half-plane packings to elements of the projective general linear 
group PGL2(Z). This is the quotient of the group GL2(Z) of invertible 2 x 2-matrices with integer 
entries, by the subgroup consisting of ±7 where I denotes the identity matrix. We also refer to the 
subgroup PSLi2(Z) consisting of those elements whose underlying matrices have determinant 1. 

Theorem 1.1. Let a, 3 > 0. There is a bijection between the set of similarities preserving the 



a b 



aa + b 



x-axis that map Vr to V a and the set of elements > £ PGL^fZ) such that = 3. The 

L c a \ " ca + a 

similarity is orientation-preserving if and only if the corresponding element is in PSLi2(Z). In 
particular, V a and Vp are similar (resp. similar via an orientation-preserving similarity) if and 

only if there exist integers a, b, c, d, with ad — be = ±1 (resp. +1), such that — - = 3. 

ca + a 

Taking a = 3 in Theorem 1 1 . 1 1 helps us to calculate the self-similarity groups. We write Symm('P) for 
the self-similarity group of the packing V, and we write Symm + ('P) for the subgroup of Symm('P) 
consisting of orientation- preserving self-similarities of V . 

Theorem 1.2. Let a > 0. Then: 

(i) IfaeQ, then V a is a strip packing and 

Symm(P a ) ^ x Z/2Z 

with subgroup 

Symm + (7 ? a) - Ax>, 
where denotes the infinite dihedral group. 

(ii) If a is quadratic over Q, then 

Symm(P a ) ^ Z. 

Let D denote the discriminant of the primitive integral polynomial with root a. Then the 
subgroup Symm + ("P Q ,) is: 

• equal to Symm("P a ) if the Pell equation x 2 — Dy 2 = —4 has no integral solution for 
(x,y); 

• the index 2 subgroup of Symm("P Q ) if x 2 — Dy 2 = —4 does have an integral solution. 

(iii) Otherwise 

Symm + (7 : ' Q ,) = Symm^a) = 1. 
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As one might imagine from the form of Theorem 1.2, there is a striking connection between half- 



plane Apollonian circle packings and continued fractions which we describe in £|4j In particular, we 
have the following result. 

Theorem 1.3. Let a,f3 > 0. The packings V a and are similar if and only if the continued 
fraction expansions of a and (3 are eventually equal (that is, become equal when initial segments, of 
possibly different lengths, are removed from each). 



As a consequence of Theorem 1.3, we obtain the following classification of the self-similar half-plane 
packings. 

Theorem 1.4. The similarity classes of self- similar half-plane (non-strip) packings are in a one- 
to-one correspondence with the finite, non-repeating sequences of positive integers, up to cyclic 
permutations. Here "non-repeating" means that the sequence cannot be realized as a concatenation 
of multiple copies of a shorter sequence. The packing has an orientation-reversing self- similarity if 
and only if the length of the corresponding sequence is odd. 



In Figures [TTJliJ at the end of the paper, we show the self-similar half-plane packings corresponding 



to the sequences (1), (2), (3) and (1,2). 

Here is a quick outline. In ^2] we give a precise definition of Apollonian circle packings and establish 
some of their basic properties. The main part of that section is then to describe a labelling system 
for the circles in a half-plane packing that are tangent to the x-axis. In §[3] we relate those labels 
to the curvatures of the circles and use this relationship to prove Theorems |1.1| and 1.2 In ^4] we 



examine the connection between half-plane packings and continued fractions which we use to prove 
Theorems 11.31 and 11.41 
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2. Apollonian Circle Packings 



We begin with a precise definition of an Apollonian circle packing. 

Definition 2.1. For the purposes of this paper, a circle in IR 2 is either a circle or a straight line, 
together with a choice of one of the components of its complement which we refer to as the interior 
of the circle. Note that what we refer to as the interior of a circle may be the unbounded component 
of its complement, and the interior of a straight line is one of the two half-planes it determines. 

An Apollonian circle packing is a collection V of circles in M? with disjoint interiors such that 

(i) there exists a set of three mutually tangent circles in V; 

(ii) if a circle C is tangent to three mutually tangent circles that are in V, then C is also in V . 

An Apollonian circle packing can be constructed recursively in the following way. 

Definition 2.2. Let be a set of three mutually tangent circles in IR 2 with disjoint interiors. 
Given pw, we define 'p( n ~\~ 1 ) to be the set of circles in M? consisting of p( n ) together with any 
circle that is tangent to three mutually tangent circles in V^ n \ The Apollonian circle packing 
generated by V^ ' is 

oo 

■p ._ j^J <p(n) 
n=0 
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It can be seen by an induction argument that the circles in "P^ have disjoint interiors, and it 



follows that V is an Apollonian circle packing in the sense of Definition 2.1 



Note that a theorem of Apollonius says that for three mutually tangent circles in M 2 with disjoint 
interiors, there are precisely two other circles tangent to all three. Each of these two circles lies in 
an interstice formed by the original three circles. See Figure [3j 

Definition 2.3. Let A, B, and C be three mutually tangent circles in M 2 with disjoint interiors. 
The complement of A U B U C in M 2 consists of five components — three of the components are 
the interiors of the respective circles, and the other two are called the interstices formed by A, B, 
and C. 




Figure 3. The two dashed circles lie in the interstices bounded by the three solid circles. 

Lemma 2.4. Let V be an Apollonian circle packing. Then V is generated, in the sense of Definition 
\2.2\ by any set of three mutually tangent circles in V . 



Corollary 2.5. IfV andV' are two Apollonian packings with a common triple of mutually tangent 
circles, then V = V . □ 



Proof of Lemma 2.4-, Since V certainly contains the packing generated by any set of three 



mutually tangent circles, it is sufficient to show that there is no room for any other circles. In 
particular, this will be true if the complement of the set of interiors of circles in V (called the 
residual set of V) has Lebesgue measure zero. A proof of this fact may be found in [U Theorem 
4.2]. □ 
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For us, the point of the recursive construction of Apollonian circle packings is that some of our 
arguments proceed by induction on the stage at which the circles are created in this process. We 
therefore make the following definition. 

Definition 2.6. Fix a generating triple V^ for the packing V. The generation of a circle C G V 
(with respect to V^), denoted by gen(C), is the unique n G Z>o such that C G ?W \ "P^" 1 ). 

We now narrow our focus to half-plane packings. Let V be a half-plane packing, that is, a packing 
that contains at least one line L. We assume that L coincides with the x-axis and that the remaining 
circles in V are in the upper-half plane. (Any half-plane packing is similar to one that satisfies this 
condition.) 

Definition 2.7. Most of our analysis of half-plane packings can be done be focusing on the circles 
in V that are tangent to the line L. We define 

V L := {C G V | C is tangent to L}. 

The 'mutually disjoint interiors' requirement of circle packings ensures that no two circles in Vl 
may be tangent to L at the same point. This property allows us to define a total ordering on the 
set Vl- We say that C is to the left of C", or C -< C, if the x-coordinate of the point of tangency 
between C and L is less than the x-coordinate of the point of tangency between C and L. In 
the case that V is a strip packing, with L' the line in V which is parallel to L, we consider the 
x-coordinate of the point of tangency between L' and L to be — oo; in other words, L' -< C for all 
C G V L with C^L'. 

Definition 2.8. Let X and Y be two tangent circles in Vl, neither of which is a line. Then 
{X, Y, L} is a triple of mutually tangent circles in M 2 and therefore determines two interstices in 
the plane. One interstice is bounded, and the other is unbounded; we refer to these as the bounded 
interstice for X and Y and the unbounded interstice for X and Y respectively. We say that 
the circle C fills the bounded (resp. unbounded) interstice for X and Y if C is the unique 
circle in the bounded (resp. unbounded) interstice for X and Y which is tangent to X, Y, and L 



(see Figure El). Note that, by Definition 2.1 C necessarily lies in V and hence also Vl- 



Remark. No circle in the bounded interstice for X and Y can be tangent to a circle in the unbounded 
interstice for X and Y . 



By Lemma 2.4 we can view the packing V as generated by the triple {X, Y, L} in the sense of 
Definition 2.2 for any pair of tangent circles X,Y £ Vl- For the remainder of this section, we fix a 
choice of X and Y and assume that X -< Y, that is, X is to the left of Y. We also assume that X 
and Y are actual circles, i.e., neither is a line. 

Definition 2.9. It is convenient to divide up the circles in Vl according to which interstice they 
are contained in. We define 

p+ = {C g V L I X < C ± Y} 



and 



V L = {C G V L I C H X or Y < C}. 



Geometrically, V^ consists of X, Y, and those circles in Vl that are in the bounded interstice for 
X and Y, while V^ consists of X, Y, and those circles in Vl that are in the unbounded interstice 
for X and Y. Note that P+U? i= V L and V£ n V£ = {X, Y}. 

Lemma 2.10. For each circle C G Vj~, C {X, Y}, there exist circles A and B in V^ of 
generation strictly less than that of C such that C fills the bounded interstice for A and B. (Recall 
that the generation of a circle in a packing V depends on a choice of generating triple; in this case, 
= {X,Y,L}.) 
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FIGURE 4. The larger dashed circle fills the unbounded interstice for X and Y, and 
the smaller dashed circle fills the bounded interstice. 

Proof. We work by induction on the generation of C. If C is generation 1, then it must be the 
circle that fills the bounded interstice between X and Y, so satisfies the lemma. Now suppose that 
gen(C) > 2. Thinking about when the circle C is added to the packing in the recursive construction 



of Definition 2.2 we see that there are exactly three mutually tangent circles of generation less than 
C that are tangent to C. One of these circl es m ust be the line L, so let A and B be the other two. 
It follows from the remark after Definition 



2.8 



that A, B G VT. Now C fills one of the interstices 
formed by A and B. We need to show that it fills the bounded interstice. 

Now exactly one of A and B must be of generation exactly one less than C. (To see this, we recall 
the procedure for recursively building an Apollonian packing. This procedure implies that as soon 
as A and B have been added, the circle C will be added in the very next generation. On the 
other hand, no two circles of the same generation are tangent since they fill different interstices.) 
Suppose this is B, so that we have gen(C) > gen(S) > 1. By the induction hypothesis, B fills the 
bounded interstice formed by two other circles of generation less than it. One of those must be A 
and let the other be D. But now we see that D fills the unbounded interstice for A and B. Since 
gen(L>) < gen(l?) < gen(C), we cannot have D = C. It follows then that C must fill the bounded 
interstice for A and B. □ 

Our main tool for keeping track of the circles in a half-plane packing V is a labelling for each circle 
in Vl by a pair of integers (a, b). The remainder of this section is devoted to the construction and 
properties of this labelling. In §3 we relate this labelling to the curvatures of the circles in V and 
use it to deduce information about similarities between different packings. 

Definition 2.11. We define a labelling function x = (x,y) : Vl — > We define the labelling 
recursively starting with x{X) = (1,0) and x{Y) = (0, 1). For the remaining circles in Vl, the label 
is determined by the following rule: 

If C is the circle that fills the bounded interstice for A and B, then 
(2.12) x(C) =x(A) +x(B). 
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At each stage of the construction of the packing from its generators, X, Y, and L, a new circle 
in Vl fihs either the bounded or unbounded interstice for a pair of circles already present. The 
equation above determines a label for each such new circle — see Figure [5j 




(a + c,b + d) 

Figure 5. The relationship satisfied by the labels. 



The main result of this section, Proposition 2.21 , tells us that the labelling function 2 is one-to-one 
and that for each pair (a, b) of coprime integers, exactly one of (a, b) and (—a, —b) is in the image 
of 2. It also gives us a necessary and sufficient condition on the labels for two circles in Vl to be 
tangent. It is convenient to start with this condition, which is stated in terms of the matrix formed 
by the labels of the two circles. 

Lemma 2.13. Let A and B be a pair of tangent circles in Vl such that A -< B. Then 



x(A) 
x(B) 



y(B) 



l. 



Proof. The proof will be by induction on gen{A,B} := max{gen(^4), gen(-B)}. The base case is 
immediate: the generation zero circles form the pair {X, Y}, which are labeled (1,0) and (0,1) 
respectively. The corresponding matrix is the identity, which has determinant 1. 

Now suppose gen{A, B} = n > 1. First, observe that we cannot have gen(A) = gen(5) = n: 
if gen(^4) = gen(U) = n > 1, then A and B were constructed to fill two disjoint interstices in 
■p( n_1 ) and cannot therefore be tangent. Hence {A,B} contains a unique circle of generation n. 
Furthermore, because a generation n circle is constructed to fill a single interstice in , p( n_1 ) ) it is 
necessarily tangent to exactly three circles of generation strictly less than n. Therefore the circle 
of generation n (either A or B) is tangent to L (generation zero), the circle in {A,B} of smaller 
generation, and a third circle C of generation strictly less than n. 

There are three possibilities for the position of C relative to A and B: C can be to the left of both, 
to the right of both, or between the two. Moreover, the generation n circle can be either A or 
B, so there are a total of six cases to consider. We only give the proof in two cases — the other 
four are nearly identical. To prove them, we use the fact that the matrix row operations of row 
addition/subtraction are determinant-preserving and that switching two rows switches the sign of 
the determinant. In each case, the final equality holds by the induction hypothesis. 

Case 1: Suppose C -< A -< B and gen(i?) = n. Then x(A) = x(C) + x(B), and 



2(A) 
x(B) 



2(A) 
2(B) - 2(A) 



2(A) 
-2(C) 



-2(C) 
2(A) 



2(C) 
2(A) 



1. 
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Case 2: Suppose A -< C ~< B and gen(A) = n. Then x(C) = x(A) + x(B), and 

1. 



x(A 
x(B 



x(A) + x{B) 
x(B) 



x(C 
x(B 



Corollary 2.14. For any circle C E Vl, gcd(x(C), y(C)) = 1. 



□ 
□ 



We now begin the proof that our labelling function x is one-to-one. We do this first for those circles 
in the bounded interstice for X and Y. 

Lemma 2.15. For C G we have x(C),y(C) > with equality only if either C = X or C = Y . 



Proof. This follows from (2.12) by induction on generation since, by Lemma 2.10, the circle C fills 
the bounded interstice of two circles of strictly smaller generation than it. □ 



Corollary 2.14 and Lemma 2.15 tell us that every circle in V~£ is labeled by a pair of nonnegative 

At the 

for Vt 



coprime integers. We now prove that every such pair is the label of a unique circle in 7^ . 
same time, we prove the converse of Lemma 



2.13 



- that if circles A,B£ have the 
determinant of the matrix formed by their labels equal to 1, then they are tangent with A -< B. 
We first need the following elementary lemma. 

Lemma 2.16. Let a and b be positive, coprime integers. Then there exist unique integers u and v 
that satisfy the following properties: 

(i) au — bv = 1, 

(ii) < u < b, and 

(iii) < v < a. 

Proof. Because a, b are coprime, we can find an integer solution {x, y) to the equation 
(2.17) ax-by = 1. 



Given a particular solution (xo,yo) to (2.17), the entire solution set is 

{(x,y) = (x + kb,y + ka) : keZ}. 

There is then a unique k G Z such that < xq + kb < b. Let u := xq + kb. Then u satisfies 
property (ii). Setting v := yo + ka, property (i) is also satisfied, and property (iii) is a consequence 
of properties (i) and (ii). 



□ 



Lemma 2.18. Let 



be a determinant 1 matrix with nonnegative integer coefficients. Then 

there exist unique circles C, C G su °h that x(C) = (a, b) and x(C) = (c, d). Moreover, C and 
C are tangent with C -< C . 



Proof. The proof is by induct ion o n max{a + b, c + d}. If max{a + b, c + d} 
. By Lemma 2.15 



1, then necessarily 

all circles in different from X and Y must have a + b > 1, 
so there can be no circles in T^i other than X and Y, labeled by the pairs (1,0) and (0, 1). 
Once we have proved that there is a unique circle with label (a, 6), we denote that circle by 



'a b 




ri 


01 


c d 







1 



At this point, therefore, we can write C, 



(1,0) 



X and C t 



(0,1) 



Y. 



Now fix an integer n > 1 and suppose we have proved the lemma, and hence constructed the circles 
C( ,6) an d Cr C)f £), for any a, b, c, d as in the 
then take a, b, c, d with max{a + b, c + d} 



C(o,6) an d C(c,<f)> for any a, b, c, d as in the statement of the lemma with max{a + b, c + d} < n. We 



n. 
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First of all, if a + b = c + d, then 

(d - b)(a + b) = d(a + b) - b(a + b) = d{a + b) - b(c + d) = ad - be = 1, 

so a + b = 1. Since c + d = a + b, this contradicts the fact that max{a + b, c + d} > 1. Therefore 
a + b > c + d or a + b < c + d. We prove the lemma in the case where a + b > c + d; the proof of 
the other case is virtually identical. 

Since a + b > c + d> 1, we have that < (a — c) + (b — d) < a + b and therefore 

max{(a — c) + (b — d), c + d} < a + b 



n. 



Since 



a — c b — d 
c d 

we show that a — c > and b — d > 0. 



1, we may apply the induction hypothesis to the matrix 



a — c 
c 



b-d 
d 



once 



Since a + b > c + d, we must have a > c or b > d. If a > c > 0, then 

ad — be = 1 cd < bc + 1 
=> cd < be 
=> d<6. 

(Note that if c = 0, then od=l, soa = d=l. It follows from the fact that a + 6 > c + d that 
6 > 1 = d.) A similar argument shows that if we assume instead that b > d, then also a > c. 

Therefore we have both a — c > and b — d > 0, so the induction hypothesis tells us that 
there are unique circles C^ a _ c ^_^ and C7 C)C rt in 7^ satisfying x (C( a -c.£>-d)) = (a — c,b — d) and 
x (C( C)( j)) = (c,d) and that, moreover, these circles are tangent with Cf a _ c ^_ d \ ~< Cr Ctd \. 

Then 



Now let C be the circle that fills the bounded interstice for C^ a _ c ^_^ and 

x(C) = (a — c,b — d) + (c, d) = (a, b), 

so there exists a circle C labeled by the pair (a,b), which by construction is to the left of and 
tangent to C( c ,d)- 

Finally, we must show that C is the only circle in that satisfies x(C) = (a, b). Suppose C" is a 
a, b). By Lemma 2.10, C" fills the bounded interstice for two circles A -< B in 

(u, u) and 

Lemma 



circle with x(C") 
V+. Then f(C") 
and B are both in P 



x(A) + x(B), so we can write x(-B) 

tells us that each of v, u, a 



l • 



< u < b and < v < a. 



2.15 



(a — v,b — u). Since ^4 
and b — u is nonnegative, so 



In fact, we have and v ^ a. Indeed, we know from Lemma 2.15 that if u = 0, then B = X. 

However, X ^ C" -< B, so X ^ B, which means we cannot have u = 0. A similar argument shows 
that a - « / 0, and so v ^ a. Therefore < u < b and < v < a, which are precisely properties 
(ii) and (hi) from Lemma 2.16 That property (i) is satisfied follows by Lemma |2. 13 since C" is 
tangent to and to the left of B. Since these three properties uniquely determine (v,u), and since 
(c, d) satisfies these three conditions by construction, we conclude that v = c and u = d. By the 



and B = C, 



(c,d) , 



following from the induction hypothesis, we have A 
and it follows therefore that C" = C. 



uniqueness of C {a _ cfi _ d) and C {c4) 



(a—c,b—d) 
□ 



Definition 2.19. Lemma 2.18 implies that any pair (a, b) of coprime nonnegative integers is the 



label of a unique circle in Vj^. As in the proof of Lemma 2.18, we denote that circle by C^ a ^y 



Lemma 



2.18 



yields a complete understanding of the labels of circles in Vt- We now define an 
operation which maps bijectively onto V~[, and use this to relate the labels of circles in to 
those of circles in Vt . 
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Let X be the unique circle which contains the three points of tangency among Cn o)> C(o,i)> an d L. 
(See Figure [6|) Define the map 



i : K 2 U {00} -»■ R 2 U {00} 



to be inversion with respect to X. 














L 



Figure 6. The inversion circle X 

Let us pause to mention some of the relevant properties of X and the map i. Note that when we 
say that i fixes a particular circle or set of circles, we mean only as sets in M 2 , not pointwise. 

(i) Inversion with respect to a circle is a bijection of order two; i.e., i o i = id. 

(ii) Since inversion maps circles to circles, i maps Apollonian packings to Apollonian packings, 
(hi) The inversion circle X intersects each of C(i 5 o)> C(o,i)' an< ^ ^ orthogonally, and therefore 



l fixes each of these three circles. Therefore, by property (ii) and Corollary 2.5 l fixes V 
(and hence Vl since L is fixed), 
(iv) The interior of X contains the bounded interstice for C(i 5 o) an d C(o,i)i an d the exterior of 
X contains the unbounded interstice for C(i,o) an d C(o,i)- Since i maps the interior of X to 
the exterior of X, and vice versa, and since property (hi) holds, it follows that t maps V 



L 



to V L and vice versa. 



Because I intersects L orthogonally, the center of X lies on L. It makes sense, then, to talk about 
a circle C £ Vl lying to the left or right of X, by which we mean that the point of tangency of C 
with L lies to the left or right of the center of X. We now record two more properties of X and t: 

(v) If C lies to the left (resp. right) of X, then i(C) also lies to the left (resp. right) of X. 
Furthermore, if C -< C both lie to the left (resp. right) of X, then l(C) ~< l(C) both lie to 
the left (resp. right) of X. 

(vi) A circle C contains the center of X (that is, the point of tangency between C and L is 
precisely the center of X) if and only if t{C) is a line parallel to L. 

As mentioned above, the reason for introducing the inversion map l is to set up a one-to-one 
correspondence between and V~[. The following lemma establishes the connection between the 
labels of circles in V~[ and their images under l, which lie in V^. 

Lemma 2.20. Let C £ VJ^ , and let = l(C) G be the image of C under the map i. Then 

x(C) = { {a '- b) V C ± C ^) 
\(-a,b) ifC m <C' 
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Proof. We will prove the statement by induction on the generation of l(C). In the case that 
gen(i(C)) = 0, we have l(C) = Cn )0 ) or t{C) = Cmu. Since Cnm and C(o,i) are fixed by t (and 
since i is one-to-one), it follows that C = C( lj0 ) or C = C( ,i)> and in both cases the statement 
holds. 

Now suppose gen(i(C)) > 1. Then, since C( a b ) = i(C) lies in V~l and is not equal to C( 10 ) or C( ,i)i 
Lemma 2.10 tells us that C( a 6 ) fills the bounded interstice for two circles C( ai)bl ), C( a2jfc2 ) 6 T 7 ^ 
with gen(C'(a 1 >6l ) ) , gen(C (a2 M ) < gen(C (a>6 )). Assume that C (aii6l ) -< C( a2) & 2 ). By definition, 
(ai,6i) + (02,^2) = (a, b). Because inversion preserves tangencies, the circles A = i{C^ ai ^) and 
B = t(C( a2 fe)) are tangent to each other as well as to C and L. 

There are a total of seven cases to consider, each corresponding to the position of the center of I 
with respect to the points of tangency of the circles C( a b ), C( ai bl ), and C( a2)b2 ) with the line L - 
the center of X could lie on one of the three circles, it could lie between two of the circles, or it 
could lie to the left or right of all three of the circles. We will prove the result for two example 
cases; the proofs in the other cases are quite similar. 

Case 1: Suppose that the center of X lies between the points of tangency of Ci ax £ x \ and Ct a ^) 
with L, as shown in Figure [7| In particular, Cr ai ^ is to the left of X and C( a {,) -< C( a2 b 2 ) are to 
the right. By property (v) above, we may conclude that A is to the left of X and B -< C are to the 
right of X. Because all three of the image circles necessarily lie in , it follows that A -< ) 
and C( !) <B -<C. 

By the induction hypothesis, we know that x(A) = (01, —61) and x(B) = (—02, 62). Since A -< B -< 
C, it follows that x(A) + x(C) = x(B), and so 

x{C) = ( - (ai + a 2 )M + 6 2 ) = (-a, b). 




Case 2: Suppose that the center of X coincides with the point of tangency between C^ a ^ and L, 
as shown in Figure [8j Then A ^ C(i,o)> C(q,i) ^ B, and C is a line parallel to L, which means that 
C is to the left of every circle in Vl- By induction, we have x(A) = (01, —61) and x(B) = (—02, 62). 
Since C -< A -< i?, we have x(A) = x[C) + x(B). It follows that 

x(C) = (01 + o a , -(61 + 62)) = (a, -6). 
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(0,6) 

Figure 8. Case 2 of Lemma l2~20l 



As mentioned above, the proofs of the other five cases are very similar to these two. 



□ 



Combining Lemmas 2.18 and 2.20 we obtain a complete understanding of how the circles in Vl are 



labeled and when two labeled circles are tangent to one another. 

Proposition 2.21. For any integers a and b with gcd(a, b) = 1, there is either a unique circle 
in Vl labelled by (a,b) or a unique circle labelled by (—a,—b), but not both. If a and b are both 
nonnegative then the label is (a, b). IfCu^s and Ci c ^\ are the unique circles in Vl with labels (a, b) 

and (c, d) respectively, then C^ a ^ is tangent to C( c ^) on the left if and only if a ^ 



1. 



Proof. Lemmas 2.18 and 2.20 



together with the fact that the inversion operation 1 acts as a 

already tells us the 'only if part 



2.13 



bijection between V~[ and V L , imply the first claim. Lemma 
of the second statement. So consider circles CV a M and Cr c d) with ad — be = 1. If all a,b,c,d are 



nonnegative then Lemma 2.18 tells us Cuj,) and CV C)C a are tangent. If one of a, b is negative, then 



the condition ad — bc = 1 implies that one of c, d must be either negative or zero. But then Lemmas 
imply that the circles fc(CV m) and t(C( C)( j)) are tangent. Since l preserves tangencies, 



2.18 



and 



2.20 



it follows that C( ,b) and C( c ^) are also tangent. That C( a ^) is to the left of C, 
Lemma 12. 131 



'M) 



follows from 
□ 



3. Self-similar half-plane packings 

We are now at a point where we may begin to describe the similarities between two half-plane 
packings. First we recall exactly what is meant by a similarity of R 2 . 

Definition 3.1. The map $ : R 2 — > M. 2 is called a similarity of M 2 if there exists some constant 
H > such that 

||$(x) -$(y)|| =n\\x-y\\ 

for all x, y G R 2 . 

Every similarity of the plane takes the form 

= fiAx + b, 
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where \i > 0, A is an orthogonal matrix, and b G M 2 . We say that $ is orientation-preserving if 
det^4 = +1, and orientation-reversing if det A = — 1. 

The set of similarities of M 2 forms a group under composition, called the similarity group of M 2 , 
which we will denote by S. The orientation-preserving similarities form a subgroup <S+. 

Similarities take circles to circles, and preserve tangency, so they take Apollonian circle packings 
to Apollonian circle packings. A key fact about the action of similarities on circle packings is the 
following. 

Lemma 3.2. Let & be a similarity of M 2 , and let A, B, and C be three mutually tangent circles 
with disjoint interiors. If A, B, and C have collinear centers, then $ is determined by the three 
circles &(A), $(B), and 5>(C), up to a reflection in the line on which the centers of the image 
circles lie. If A, B, and C have non-collinear centers, then is completely determined by $(A), 
<S>(B), and $(C). 

Remark. Here we mean that the 'center' of a line L lies 'at infinity' orthogonal to L in the direction 
of its chosen interior. If one of the circles A, B, C, say A, is a line, then the collinearity condition 
is satisfied if and only if another of the circles, say B, is also a line, parallel to A. In this case $ is 
determined up to a reflection in the line through the center of C that is orthogonal to A and B. 

Proof. Since similarities form a group, it suffices to consider the similarities that fix A, B, and C 
(as sets, not pointwise). If a similarity fixes the circles A, B and C, then it fixes their centers. A 
similarity that fixes three non-collinear points must be the identity. One that fixes three distinct 
collinear points is either the identity or a reflection in the line formed by them. □ 

Definition 3.3. For $ £ 5 and an Apollonian packing V, we write 

$.p:={$(C) :CeV}. 

Two packings V and V are similar if V' = $ ■ V for some $ 6 S. 

The group 

Symm(-P) := {$ € S | $ • V = V} 
is the self-similarity group of V . The subgroup of Symm('P) consisting only of orientation- 
preserving similarities is the orientation-preserving self-similarity group of V, denoted by 
Symm + ("P). A packing V is self-similar if Symm('P) is nontrivial. 

In order to establish similarity between two packings, we look at the curvatures of the circles 
involved. 

Definition 3.4. For a circle C in M 2 , the curvature of C, denoted curv(C) is the reciprocal of 
the radius of C. A straight line in R 2 is considered to have curvature zero. 



Lemma 3.2 allows us to check similarity by looking only at the curvatures in a triple of mutually 
tangent circles in each packing. 

Lemma 3.5. The packings V and V are similar if and only if they contain triples of mutually 
tangent circles (A,B,C) and (A' , B' ,C) respectively, such that there exists \i > with 

curv(^4') = /icurv(A), curv(S') = ^curv(S), curv(C') = //curv(C). 

Proof. If V and V 1 are similar via similarity $ with scale factor /i, then take any triple (^4, B, C) 
and set A' = $>(A),B' = &(B),C' = $(C). To prove the converse, choose a similarity $ of the 
plane that takes A to A', B to B' and C to C . (One can choose a translation composed with 
dilation to get A to A' , add a rotation to get B to B' , then add a reflection if necessary to get C to 
C.) By construction, V' and $(V) both contain the triple {A',B',C}, and therefore V' = $(V) 



by Corollary 2.5 □ 
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Turning now to half-plane packings, that is, those that have a straight line for at least one of the 
circles, recall that we can focus on the following packings. 

Definition 3.6. Let a G R. Then V a is the packing generated by a triple {X, Y, L} of mutually 
tangent circles, where L is the x-axis in M 2 , X is a circle of curvature a 2 tangent to (and above) L 
at the origin, and Y is a circle of curvature 1 resting on L and tangent (on the right) to X. The 
generating triple is illustrated in Figure [2] in the Introduction. 

Lemma 3.7. Every half-plane packing V is similar, via an orientation -preserving similarity, to T^q 
for some a > 0. 



Proof. Choose any two tangent circles in V that are tangent to a line but are not themselves lines. 
Taking /x to be the ratio of their curvatures (in the appropriate order), and a = y/Ji, this follows 
from 13.51 □ 



The key to analyzing Apollonian circle packings is the following result, due to Descartes. This 
describes the relationship between the curvatures of four mutually tangent circles in the plane. A 
selection of proofs of this are given in [9]. 

Theorem 3.8 (Descartes' Circle Theorem). Let w, x, y, and z represent the curvatures of four 
mutually tangent circles in the Euclidean plane. Then 

2(w 2 + x 2 + y 2 + z 2 ) = (w + x + y + z) 2 . 

For half-plane packings, we apply this Theorem in the case where one of the four circles is a line, 
i.e. has zero curvature. In this case, the quadratic relationship boils down to a linear relationship 
between the square roots of the curvatures of the circles. 

Corollary 3.9. Let a 2 , f3 2 and r y 2 represent the curvatures of three mutually tangent circles all 
tangent to a line L, where a > (3 > and 7 > 0. Then 

(3.10) j = a±/3. 

In particular, when the circle of curvature 7 2 lies in the bounded interstice formed by the others, 
we have 

(3.11) 7 = a + /3. 



Proof. The proof of (3.10) follows from Theorem 3.8 by setting w = and applying the quadratic 



formula appropriately. The proof of (3.11) follows from the fact that the circle in the bounded 



interstice has a curvature at least as large as that of the two circles surrounding it. □ 



We illustrate (3.11) in Figure 



Observe that the illustration of this equation is virtually identical 
which shows the recursive labelling process defined in ^2] It is this 

A key 



to the illustration in Figure p5 

linear relationship between the curvatures of tangent circles that inspires that labelling. 



consequence of this connection is Lemma 3.13 below. 



Definition 3.12. Consider now some fixed real number a > and recall the packing V a from 
Definition 3.6 We write V a ,L f° r the set of circles in V a that are tangent to the line L (that is, 
the x-axis). Let C(i,o) denote the circle of curvature a 2 that is tangent to L at the origin, and let 
C(o,l) denote the circle of curvature 1 that is tangent to L and to C( 10 ) on the right. As described 
in the previous section, these choices determine a unique label for each circle in V a ,L- When we 
need to specify the underlying a we use a superscript, as in C" ab y but we often drop the a when 
context allows. 
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a + P 



Figure 9. Descartes' Theorem for circles in Vl- The circles are labelled by the 
square roots of their curvatures. 



Lemma 3.13. The circle C^ a ^ in V a ,L h> as curvature given by 

curv(C( aj6) ) = (aa + b) 2 . 

Moreover, aa + b > 0. 

Proof. The proof is by induction on the generation of C( a m with respect to the generating triple 
{L, C(i,o)j C(o,i)}- The result is immediately seen to hold for the generation zero circles Cno) and 
CV 0i i), since they were chosen to satisfy curv(C( 10 )) = a 2 and curv(C( 0j i)) = 1. 

Now suppose gen(C( a b)) = n > 1. The circle C( a ,fe) was constructed to fill an interstice bounded by 
three circles of generation strictly less than n; since C^ a ^ is tangent to L, L is necessarily one of 
those circles. Because the other two circles are tangent to L as well, we can call them C( ai 5^ and 



C(a 2 ,b 2 ) with C (ai)6l ) -< C (a2i62 ). 



We do the case where C( a ^ fills the bounded interstice between 



C(ai,bi) an d C(a 2 ,fe 2 )- The case where it fills the unbounded interstice, either to the left or right, is 
similar. 



By Definition 2.11 we have a = a\ + 02 and b = bi + bi- By Corollary |3.9| then, we get 



Y / 'curv(C (a>6) ) = ^/curv(C (ai)6l) ) + ^curv(C (o2ib2 )) 
= aia + 61 + 02a; + &2 
= aa + b, 

where the second equality holds by induction. 



□ 



The following is an immediate consequence of Lemma 3.13 

Corollary 3.14. If a $ Q, then no two circles in V a have the same curvature. 



Proof. Suppose curv(C( a6 )) = curv(C( a / b /)). Then, by Lemma 
Since a is not rational, the only way for this equation to hold is for a 
uniqueness statement in Proposition 2.21, implies that C( aj , 



3.131 we have aa + b 
a' and b 



a' a + b'. 
b' which, by the 
□ 



We are now in a position to prove our first main result, identifying the set of similarities between 
the two packings V a and Vp when a, j3 are positive real numbers. This is Theorem |1.1| from the 
Introduction. We start by showing how to associate a matrix to such a similarity. 

Definition 3.15. Fix a, ft > and let $ be a similarity of M 2 such that $ • Vp = V a . Also assume 
that ^(L) = L, where L is the x-axis, that is, the chosen line in each packing. Then <3? takes 
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mutually tangent circles in V@,l to mutually tangent circles in V a L- In particular, we have 



c; 



(0,6) > 



(0,1)' 



for some integers a,b,c,d. By Lemma |2. 13 



a 6 
c d 

is an integer matrix of determinant ±1. The determinant is +1 if C^ ab ^ -< C(£d)> i n which case $ 
is orientation-preserving, and — 1 if <& is orientation-reversing. We denote this matrix by A($). 

Theorem 3.16. Let a,B > 0. The construction A of Definition 3.1^ determines a bijection 
between the set of similarities 0/M 2 that take Vp to V a (and fix the x-axis) and the set of matrices 
a b~ 



d 



in PGIj2(Z) that satisfy 



aa + b 



ca + d 

Furthermore, the restriction of A to orientation-preserving similarities is a bijection onto the set 
of elements o/PSL2(Z) = SL2(Z)/{±1} with this property. 

Proof. We first show that A(<&) satisfies the condition that B = Because $ is a similarity, 



ca 



d 



there exists some A > such that curv($(C)) = Acurv(C) for all C G Vp. Since curv(C^ Q ^ 



B 2 



and curv(C^ 01 j) = 1, it follows that curv(C^^) = A/3 2 and curv(C^^) = A. By taking square 
roots and applying Lemma 3.13 we may conclude that 

aa + b = vA/3 

ca + d = VA, 

aa + b 



which we may rewrite as 



ca + d 



To show that A is injective, suppose A($) = A($') in PGL 2 (Z). Then, if $(Cf 1Q) ) = C^ b) and 
^'(^(l 0)) = C( a > (,')> we must have (a, b) = ±(a', V). But, by Proposition 2.21 , only one of (a, b) and 



(—a, —b) is the label of a circle in V a . Therefore, in fact (a, b) = (a', V) and so $(C^ ^) = 3>'(C£ Q s). 



Similarly x) ) = $'(cf 



(0,1)' 



Since also = <&'(L), Lemma 



3.2 



tells us that $ = 



Now let 



a b 
c d 



be an element of PGL2(Z) such that B = — -. Because the determinant of this 

ca + d 



matrix is ±1, Proposition 2.21 tells us that either (a, b) or (—a,—b), but not both, is the label 
of a circle in V a ,Lj and that the same holds for (±c, ±d). Furthermore, these circles are tangent. 
Because we can multiply the matrix by —1 and not change it in PGL2(Z), we may assume that 
(a, b) is the label for a circle in V a L- Now we need to show that (c, d) is also the label of a circle 
in VaL- Suppose that (— c, — d), rather than (c, d), is a label in V a L- Lemma 3.13 tells us that 

aa + b > 0. This statement, along with the fact that — - = B > 0, implies that ca + d > 0, 



ca 



and therefore (— c)a + (— d) < 0, which contradicts Lemma 



may therefore conclude that C? b ^ and C? c d ^ form a pair of tangent circles in V a ,L- Then, since the 



d 



3.13 



applied to the circle C 



(%,-*>■ We 



ratio of the curvatures of C, 



(0,6) 



and C? _n is 



M) 



(aa + b) 2 
{ca + d) 2 



/? 2 
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which is the same as the ratio of the curvatures of C; 
similarity $ between V a and Vg such that 

Km = 



(1.0) 



and C: 



P 

(0,1) ■ 



Lemma 



3.5 



determines a 



This shows that A is a bijection. We have already noted that A(<£) G PSL 2 (Z) if and only if $ is 
orientation-preserving, which gives us the last part of the Theorem. □ 

We have the following corollaries. 

Corollary 3.17. V a is a strip packing if and only if a G . 



Proof. Fi rst, w e note that V a is a strip packing if and only if it is similar to the packing V\. By 
this is true if and only if there is an integer matrix ^" ' 



Theorem 
a- 1 

that 



3.16 



d 



of determinant ±1 such 



a + b 



c • 1 + d c + d 

2 with p, q > and gcd(p, q) 



suppose a 
aq — cp = 1, and set b 



a. Certainly, if such a matrix exists, then a is rational. Conversely, 
1. Let a and c be positive integers that satisfy 

P 



p — a, d = q — c. By construction, 



a b 
c d 



1 and 



a + b 



c + d q 



a. 



□ 



Corollary 3.18. If V a is self-similar, then a is the root of a quadratic polynomial with rational 
coefficients. 



3.16 

that 



this corresponds to a nontrivial element 



Proof. If V n i s self-similar, then there is a nontrivial similarity <3? that maps V a to itself. By 
Theorem 
i.e., such 

ca 2 + (d — a)a — 6 = 0. 



G PGL 2 (Z) such that 



aa + b 
ca + d 



It is easy to check that the only way for all three coefficients to be zero is for a = d = ±1, b 

'a b~ 



which contradicts the fact that the matrix 



d 



is not the identity in PGL2(Z). 



0, 

□ 



The rest of this section is concerned with proving the converse of Corollary |3. 18 if a is the root of a 
quadratic polynomial with rational coefficients, then V a is self-similar. This follows from Theorem 



3.22 below, which is Theorem 1.2 of the Introduction 



Theorem 3.22| goes beyond identifying which packings are self-similar. We in fact calculate the 
self-similarity groups of all the packings. To do this we first show that, in the case (3 = a, the 



bijection of Theorem 3.16 is a group isomorphism. This is the content of the following proposition. 



Proposition 3.19. Let a > be irrational. Then there are group isomorphisms 

aa + b 



Symm("P Q ) = Stab(a) := 
Symm + ('P Q ) = Stab+(a) := 



a b 
c d 

a b 
c d 



€ PGL 2 (Z) 
G PSL 2 (Z) 



ca + d 
aa + b 
ca + d 



a > , and 



a 



Proof. Since a is not rational by assumption, Corollary |3.17 tell us that V a is not a strip packing, 
so the x-axis L is the unique line in the packing V a . Therefore every self-similarity of V a maps L 
to L. If we take (3 = a, then Theorem 3.16| states precisely that we have bijections of the above 
forms given by the construction A. Now we show that when f3 = a, these bijections are group 
isomorphisms. 
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A($) 



a b 
c d 



and A($') 



6' 



3.14 



and 



We must show that A(d> o **) = [« °j [5 = [£ Tdc' S TSJ " B * C ° ro11 ^ 
Lemma |3.2[ it will suffice to show 

(3.20) curv(($ o $')(C (lj0 ))) = ((aa + bc')a + (a& 7 + W)) 2 and 

(3.21) curv(($ o *')( c (o,i))) = ((<»' + dc> + (cb' + tM')) 2 - 

In fact, it will suffice to show only that ( |3.21 ) holds: since <!> o is a similarity, we must have 
curv((<£ o $0(C(i,o))) = Acurv (C(i,o)) = Aa^ and curv((<I> o $')(C(o,i))) = A eurv(C( 0)1 )) = A for 
some A > 0; therefore, if (3.21) is satisfied, the fact that the product of the matrices is still an 
element of Stab(a) will force ( |3.20| ). 

The scale factor A under the composition <3? o $' is the product of the scale factors /i and // under 
the maps and respectively. Since curv(C(o,i)) = 1, it follows that 

ji = curv($(C( 0) i))) = curv(C( c = (ca + d) 2 and 

m' = curv ( $/ (C(o,i))) = curv(C (c / )d /)) = (ca + d') 2 . 

Therefore we will have shown that A is a group homomorphism if we can show that 

(ca + d)(c'a + d!) = (ca' + dc')a + (cb' + dd'). 

Indeed, 



a' a + b' 
da + df 



a 



ca 2 + (d' — a')a — b' = 

cc'a 2 + (cd! — ca')a — cb' = 

cc'a 2 + (cd' + dc — ca! — dc')a — cb' = 

cc'a 2 + (cd' + dc')a = (ca + dc')a + cb' 

cc'a 2 + (cd' + dc')a + dd' = (ca' + dc')a + (cb' + dd!) 

(ca + d)(c a + d') = (ca + dc')a + (cb' + dd'). 



□ 



Now we can prove the main result of this section. 
Theorem 3.22. Let a > 0. Then: 

(i) If a € Q, tfiera 7> a is a strip packing and 

Symm(-P Q ) ^D m x Z/2Z 

ure£/i subgroup 

Symm + (P Q ) ^ ZA^,, 
where denotes the infinite dihedral group. 

(ii) If a is quadratic over Q, then 

Symm(P a ) Z. 

Let D denote the discriminant of the primitive integral polynomial with root a. Then the 
subgroup Symm. + (V a ) is: 

• equal to Symm(7 : ' a ) if the Pell equation x 2 — Dy 2 = —4 has no integral solution for 
(x,y); 
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• the index 2 subgroup of Symm('P Q ) if x 2 — Dy' L 
(iii) Otherwise 

Symm + (V a ) = Symm('P Q ,) = 1. 



4 does have an integral solution. 



Proof. First of all, it follows immediately from Corollary |3.18| that, for any a > that is neither 
rational nor quadratic, V a is not self-similar. Therefore Symm('P a ) = Symm + ("Pa- ) = 1 for all such 
a. 

Now suppose that a is rational. Then Corollary |3.17 tells us that V a is a strip packing. In this 
case, the full self-similarity group is generated by a translation 'along' the strip, a reflection in 
a line perpendicular to the strip, and a reflection that interchanges the two lines. The resulting 
group is isomorphic to x Z/2Z where is the infinite dihedral group. The subgroup of 
orientation-preserving self-similarities of V a is generated by the translation and the rotation given 
by combining the two reflections. This subgroup is isomorphic to -Doo- 

The main focus of our work is the case where a is of degree precisely 2 over Q. We have already 
shown in Proposition 3.19 that 

Symm(P a ) 5 



a b 
c d 



a b 
c d 



E PGL 2 (Z) 
G PSL 2 (Z) 



aa + b 
ca + d 

aa + b 



a 



ca + d 



a 



and that 

Symm + (P a ) 5 
so it suffices to calculate these stabilizer groups. 

It turns out that the elements of PGL2(Z) that fix a are closely related to the solutions to the Pell 
equations 

(3.23) x 2 - Dy 2 = ±4, 

where D = q 2 — Apr is the discriminant of the primitive integer polynomial f(x) 
satisfied by a, with p > 0. 



px + qx + r 



Define the sets Q and Q+ as follows: 

x + yy/D 



Q :-- 



and 



x + y\/~D 



x,y 6 Z, x 2 -Dy 2 = ±4 



x, y 6 Z, x — Dy 



where D is as in the previous paragraph. One can easily check that Q is a group under multiplication 
with subgroup Q+. Furthermore, one can show (see [8l Theorem 1.9], for example) that 

(3.24) 5 = Zx{±l}. 



Our proof of Theorem 3.22 is given by relating the stabilizer of a in PGL^Z) with the group Q. 

Recall that the integers p, q, r are the coefficients of the primitive integer polynomial satisfied by a 
and that D = q 2 — Apr. We first construct a group homomorphism 

aa + b 



r-.g 

x + yy/~D 



a b' 
c d 

x-yg 

2 



E PGL 2 (Z) 



ca + d 



n 



VP 



-yr 

x+yg 

2 J 
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The matrices in the image of T consist of integer entries because: 

x ± yq = x 2 — y 2 q 2 (mod 2) 

= x 2 — y 2 (q 2 — 4pr) (mod 2) 



x 2 -Dy 2 = ±4: = (mod 2). 



These matrices stabilize a because: 
pa + qa + r = 



ypa + yqa + yr = 



ypa + 



x + yq x — yq 



a + yr = 



x—yn 

2 yH a — yr 
ypa + ^ 



a, 



and are invertible because 



x —yq _„. r 

2 y 

yp ^ 



x 2 -y 2 q 2 2 
Yy pr 



1 

4 
±1 



x 2 - Dy 2 ) 



where the sign is positive if and only if (x, y) satisfies x 2 — Dy 2 = +4. To see that T is a group 
homomorphism, we check: 



'x + yy/D x' + y'VD^ 



xx'+Dyy' , xy'+x'y 
n r o 



xx' +yy' (q —4pr)—xy' q—x'yq —xy'r—x'yr+yy'qr—yy'qr 



x' yp+xy' p+yy' pq—yy' pq xx 1 '+yy' ' (q 2 ' —ipr)+xy' ' q+x' 'yq 
2 4 . 



x—yn 

-yr 

yp 



2 

y'p 



-y'r 



r / x + yVD \ r / x' + i/y/D 



Combining the map T with the isomorphism of Proposition 3.19 we have now shown how to con- 



struct, for each solution to (3.23), a self-similarity of V a . To prove our Theorem, we calculate 



the kernel and image of the homomorphism T. First, we show that F is surjective, which implies 



that every self-similarity of V a arises from a solution to (|3.23|) in the manner described above. So 

"a b~ 



suppose we are given a matrix A 



c d 



that stabilizes a. In particular, it follows that 



ca 2 + (d — a) a — 6 = 



This polynomial is therefore an integer multiple of the primitive polynomial px 2 + qx + r with root 
a. That is, there exists m G Z such that 



c = mp 
d — a = mq 
—b = mr. 
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Now set 

x = a + d 
y = to. 

We clearly have x, y G TL and 

x 2 - Dy 2 = (a + d) 2 - {q 2 - Apr)m 2 
= (o + d) 2 - (a - d) 2 - Abe 
= A{ad - be) = ±4. 

so x+ y^ £ g. (Moreover, this is in Q + if and only if A G PSL^Z).) It is easy to check that 
T{ x+y ^~^ ) = A as required. 

Finally, a:+ ^ v/ ^ i s i n the kernel of V if and only if 

x — yq = x + yq = ±2 , yr = yp = 0. 
Since p cannot be zero (a ^ Q), it follows that y = 0, and therefore x = ±2. In other words, 

ker(r) = {±1}. 

Putting together our various isomorphisms and using ( |3.24 ), we now have 



Symm(P a ) £/{±l} ^ Z. 

We have also seen that the orientation-preserving self-similarities correspond under this isomor- 
phism to the subgroup C/+/{±1}. There are two possibilities here. One is that the generator for 
Q is in Q + . In this case the groups are equal and all the self-similarities of V a are orientation- 
preserving. This happens when there are no integer solutions to the equation 

x 2 - Dy 2 = -A. 

The other possibility is that the generator z = x ° + yovD £ Qr g - g no ^. - n g u ^. however, z 2 

is in Q + and so Q + is an index 2 subgroup of Q. In this case, Sym.m. + {V a ) is an index 2 subgroup 
in Symm("P Q ) as claimed. □ 

Corollary 3.25. The half-plane packing V a is self-similar if and only if a is rational or quadratic 
over Q. In the quadratic case, V a is self -similar via an orientation-reversing self- similarity if and 
only if the equation x 2 — Dy 2 = —A has an integral solution (x, y), where D is the discriminant of 
the primitive integral polynomial with root a. 

A. Half-plane Packings and Continued Fractions 

In this section, our main goal is to describe how the continued fraction of a positive real number 
a manifests itself geometrically in the half-plane packing V a - Recall that the standard continued 
fraction expansion of a positive real number a is a representation of the form 

1 

a = ao H 



1 

a\ -I 



1 

a 2 H 

for some integers a^, with oq > and > for all k > 0. It is also standard (and more practical 
typographically) to express this expansion simply as a = [do, a±, . . .]. 



APOLLONIAN CIRCLE PACKINGS OF THE HALF-PLANE 



23 



To begin, we recall the algorithm for computing the continued fraction expansion of a positive real 
number a. The continued fraction expansion is computed by successive iterations of the following 
algorithm, which we refer to as the continued fraction algorithm. The input for the algorithm 
is the number ao = a. Each step of the algorithm takes a n and determines a n +i- 

(A) If a n > 1, let a n+ i = a n - 1. 

(B) If < a n < 1, let a n +i = 

(C) If a n = 0, halt. 

Recording the sequence of steps obtained when applying this algorithm to a positive real number 
a we get something like 

ABAABAAC. 

The positive integer a>k from the continued fraction expansion corresponds precisely to the length of 
the (k + l)th string of consecutive A's. For example, the above sequence represents the application 
of the continued fraction algorithm to «o = f • The resulting continued fraction expansion is 

[1,2,2] = ! + ^=^. 
2+ 2 

The sequence (a n ) in this case is: 

(«) 7(^2(^5(^3(^1^ <^ 

v ; 5 5 2 2 2 

Notice that the continued fraction expansion for a n is the same as that for a, but with an 'initial 
segment' removed. For example, if «o = [2, 3, 4, 5, 6, 7], then we have 

ai = [1,3,4,5,6,7] 
a 2 = [0,3,4,5,6,7] 

a 3 = [3,4,5,6,7] 

a 4 = [2,4,5,6,7] 

a 5 = [1,4,5,6,7] 

a 6 = [0,4,5,6,7] 
a 7 = [4,5,6,7] 

and so on. 

Turning now back to Apollonian circle packings, we define a circle replacement algorithm. 
The input of this algorithm is an ordered pair (Xq,Yq), where Xq and Y$ are tangent circles in a 
half-plane packing V that are also tangent to a chosen line L £ V . We also require that Yq is not 
itself a line. At the (n + l)th step of the algorithm, we replace the pair (X n ,Y n ) with a new pair 
of circles (X n+ i,Y n+ i): 

(A) If curv(X n ) > curv(y n ), take Y n+ \ = Y n and take X n+ \ to be the circle that fills the 
unbounded interstice for X n and Y n (in the sense of Definition 2.8). Note that Corollary 



3.9 implies that 

A/curv(X n+ i) = A/curv(X n ) - 1/ curv(Y n 

(B) If < curv(X n ) < curv(l^), take X n+ \ = Y n and Y n+ \ = X„ 

(C) If curv(X n ) = 0, halt. 
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As with the continued fraction algorithm, we are interested in the sequence of steps involved when 
the algorithm is performed to a given starting pair of circles. (For example, we might obtain the 
sequence AAB AB AAAAB AAC '.) Our main observation is then the following. 

Lemma 4.2. Let a be a positive real number, and let {Xq, Yq) be the two circles used to construct the 
half-plane packing V a : Xq and Yq are tangent to each other and to the x-axis L, and curv(Ao) = a 2 , 
curv(Yb) = 1- Then the sequence of steps (A, B, or C) performed in applying the continued fraction 
algorithm to a is the same as the sequence of steps performed in applying the circle replacement 
algorithm to (Xq,Yq). Moreover, we have 

1 /curv(A" n ) 



a. 



y/curv(Y n ) 



for all n > 0. 



Proof. The proof is by induction on n. For n = 0, this is the claim 

v / curv(X ) 



a 



v /curv(y ) 
which is true by the choice of Xq and Yq. 

Suppose that the claim holds for a n and (X n , Y n ). Then a n > 1 if and only if curv(X n ) > curv(Y n ) 
and a n = if and only if curv(X n ) = 0. This tells us that the next step (A, B, or C) will be the 
same for both algorithms. So it remains only to verify that the formula still holds for a n +i and 
{X n+ \, Y n+ \). 

Suppose that a n > 1. Then we have a n +i = a n — 1, so it is sufficient to show that 

v / curv(A ra+ i) _ v / curv(X n ) _ 
v /curv(y n+ i) 1 /curv(l" n ) 
We have Y n+ \ = Y n , so it is enough to show that 

i/curv(A" n+ i) = ■ v /curv(X n ) - y/cvrv(Y n ) 
which follows from Corollary |3.9| as mentioned above. 
Finally, suppose that < a n < 1. Then 

1 _ v / curv(y n ) _ v / curv(X n+ i) 



a n .fi 

ct n y / curv(X n ) A /curv(y n+ i) 



□ 



Figure 10 shows the circle replacement algorithm applied to the packing V a for a = I. The circles 



are labeled by the square roots of their curvatures. Compare this to (4.1) as an illustration of 
Lemma 14.21 



Recall that the continued fraction expansion of a real number a determines a sequence of rational 
numbers ^ that converge to a. These are the convergents of a and are given by truncating the 
continued fraction expansion of a. Thus if 

a = [00,01,02,...] 

then set 

Pn r , 

— :— ao, 01, . . . , a n , 

where p n and q n are nonnegative coprime integers. They satisfy the recurrence equations 

(4.3) p n = p n -2 + CLnPn-1, q n = + a n q n -\. 
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Figure 10. The circle replacement algorithm applied to the packing Vr ■ The circles 

5 

are labeled by the square roots of their curvatures. 

We now observe that the convergents of a appear in the labels (in the sense of Q of the circles in 
the circle replacement algorithm applied to the packing V a . 

Lemma 4.4. Let a be a positive real number. Let (Xq,Yq) = (Cmq), C(o,l)) ^ e the generating 
circles for the packing V a . The sequence of distinct circles in the sequence (Yj) defined by the 
circle replacement algorithm is 

u (0,l); ^(qo,-po)-> U (-<?i,Pi)' °(<?2,-P2)> U (-<?3,P3)' 

In particular, if a Q, then 

lim curv(Yj) = 0. 

j'->°o 

Proof. We have Yq = C( ,i) an d Xq = Cn y The first new Yj will appear after the first application 
of step (B) of the algorithm, that is, after ao + 1 steps. At this point we have 



Ya +1 — X aQ — C(! _ aQ ) — C( go _ pQ ) 



and 



X ao +i — Y ao — C( 0) l). 

Now suppose, inductively, that immediately after the nth application of step (B) we have 
Running the algorithm until after the next application of (B), that is a n + 1 times, we have 

^AT+a„+l — ^iV+a n — C(-q n „2-a.n<Jii-l,Pn-2+cinPn-i) — ^(—q n ,p„) 

and 

The first claim now follows by induction on n. 



By Lemma 3.13 we have 



cur 



v(C( ±gniT p n) ) = (q n a -p n ) =q n [a 



Pn 
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A basic fact about the convergents for continued fractions [7J Theorem 171] is that 

Pn 



a 



1 

< — 



so 



1 

curv(C (±(?niTPn) ) < -2 . 



In 



It follows from (4.3) that q n — > oo as n — > oo, so the curvatures of the Yj tend to zero. □ 



Our goal is now to use this relationship between the circle replacement and continued fraction 
algorithms to give new criteria for two packings to be similar, and a new way to understand the 
self-similarities of a given packing, both in terms of continued fraction expansions. To do this we 
have to know that the circle replacement algorithm involves 'enough' of the circles in the packing 
to be able to detect any similarity between two packings. The following lemma is key to this. 

Lemma 4.5. Let V be a half-plane packing (but not a strip packing). Fix an ordered pair of circles 
(Xo, Yq) as in the definition of the circle replacement algorithm. Let X and Y be any pair of tangent 
circles in Vl such that Xq and Yq are contained in the bounded interstice formed by X and Y . Then 
one of the pairs (X, Y) and (Y,X) appears as (X n ,Y n ) in the application of the circle replacement 
algorithm to (Xq,Yq). 

Proof. The second part of Lemma |4.4| implies that the radii of the circles X n and Y n increase 
without bound as n tends to infinity. Therefore there is some smallest integer N such that Xn 
does not lie in the bounded interstice formed by X and Y. Since Xn is tangent to Xjv-i, the 



remark following Definition 2.8 implies that Xn also cannot lie in the unbounded interstice for X 
and Y. Hence Xn is equal to one of X,Y, say Y without loss of generality. 

Now Yjv is either equal to X n for some n < N, or is equal to Yq. Either way, Yn is in the bounded 
interstice formed by X and Y and so, in particular, is smaller than Xn- This tells us that the 
next step in the circle replacement algorithm is (B); i.e., Xn+i = Yjy and Ijv+i = Xn = Y. Since 
Yn+i is larger than Xn+i, we next repeat step (A) until Xn+k is larger than Yn+k = Yn+\ for 
some K > 0. Then Xn+k is not in the bounded interstice formed by X and Y, but Xn+i is. This 
means we can find a smallest M with iV + 1 < M < N + K such that Xm is not in the bounded 
interstice formed by X and Y. Since it is tangent to Xm—i, this circle Xm also cannot be in the 
unbounded interstice, so must be one of X and Y. But it is not Y since Ym is. Therefore we have 
Xm = X and Ym = Y which completes the proof. □ 

We can now relate properties of the continued fraction expansion of a positive real number a to 
geometric properties of the half-plane circle packing V a - 



4.1. Strip packings. We already saw in Corollary 3.17 that V a is the strip packing if and only if 



a € Q. This is now reflected in the fact that the continued fraction expansion for a halts if and 



only if a € Q. We can see from Lemma 4.2 that the continued fraction expansion of a halts exactly 



when the corresponding circle replacement algorithm produces a circle of curvature 0, that is, a 

7 
5 



straight line. This is illustrated in the example of a = I displayed above. 



4.2. Similar packings. We can determine whether the packings V a and Vp are similar by exam- 
ining the tails of the continued fractions of a and (5. 

Definition 4.6. Let us say that a and /3 have eventually equal continued fraction expansions 
if there is some k,N £ Z such that a n = 6 ra +fc for all n > iV (where [at] is the continued fraction 
expansion of a and [bj] is the continued fraction expansion of j3). For example, [1, 2, 3, 4, 4, 4, 4, ... ] 
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and [2, 7, 4, 4, 4, 4, ... ] are eventually equal. For the purposes of Theorem 4.7, we say that two finite 



continued fraction expansions are eventually equal since both expansions terminate. 

It is easy to see that a and (3 have eventually equal continued fraction expansions if and only if 
there exist m, n > such that a n = (3 m (where these are the sequences obtained by applying the 
continued fraction algorithm to a and f3). This observation yields the following classification of 
half-plane packings up to similarity. 

Theorem 4.7. For positive real numbers a, j3, the circle packings V a and Vp ore similar if and 
only if a and f3 have eventually equal continued fraction expansions. 

Proof. A quick proof of this result follows by identifying each of the conditions in the statement 
with the condition that there exist p,q,r,s G Z with ps — qr = ±1 and ^a+g = /3. For the circle 



packings, this is Lemma 3.16. for the continued fractions, it is fTJ Theorem 175]. However, we give 
a more interesting proof arising from the direct comparison between the continued fraction and 
circle replacement algorithms. 

First, note that a and (3 have finite continued fraction expansions if and only if a and /3 are rational, 



which is equivalent by Corollary 3.17 to V a and Vp both being strip packings, which are similar. 
We may therefore assume that a and (3 are irrational. 

Suppose a and /3 have eventually equal continued fraction expansions. Then a n = f3 m for some 
m,n. This means that the ratio of the curvatures of a pair of tangent circles in V a , both tangent 
to L, is equal to the ratio of the curvatures of a pair of tangent circles in Vp, both tangent to L. It 



follows by Lemma 3.5 that there is a similarity between V a and Vp. 

To prove the converse, suppose V a and Vp are similar. Then there is a pair of circles (X' q ,Yq) 
in V a whose ratio of curvatures is equal to (3 2 , in addition to the original pair of circles (Xq, Yq) 
in V a whose ratio of curvatures is a 2 . The key step is the following claim: if we apply the circle 
replacement algorithm to each of these pairs of circles, they will eventually coincide; that is, there is 
some pair of circles (X, Y) in V a that appears both as (Xm, Ym) and [X' N , Y^) for some M, N e N. 
Note that the circle replacement algorithm only sees ratios of curvatures and not the curvatures 
themselves, so the circle replacement algorithm will generate the same numerical data for (X' , Y') 
as it would for the corresponding pair of circles in the packing Vp. From this claim, it follows that 
a and /3 have eventually equal continued fraction expansions since 

v / curv(X M ) y/curv(X) ^/curv{X 1 N ) 
v / curv(y M ) Vcurv(Y) V c urv(Y^) 

Now let us prove that the circle replacement algorithms corresponding to a and /3 eventually 
coincide in the sense described in the previous paragraph. Let (X n ,Y n ) be the pairs obtained from 
applying the circle replacement algorithm to (Xq,Yq), as defined in the previous paragraph. We 
suppose without loss of generality that Xq -< Yq and X' -< Y '. If one of the pairs contained the 



other in its bounded interstice, say Xq ^ X -< Y ' ^ Yq, then by Lemma 4.5 there would be M 
such that {Xm,Ym} = {X ,Yq}. Suppose instead that 

Xq -< Yq ± X'q -< Yq 1 . 

Suppose that there is no pair (X n ,Y n ) that contains (X ,Yq) in its bounded interstice. Then, 
for each n, one of the circles X n and Y n , has its point of tangency with the x-axis between the 



corresponding tangency points of Yq and X' . By the second part of Lemma 4.4 this means that 
there are arbitrarily large circles, all disjoint, with tangency points in this fixed interval. A little 
geometry shows that if two disjoint circles of radii R and R' are tangent to the x-axis, then their 
points of tangency are at least 2\/ RR' apart. This gives us a contradiction and so we deduce that 
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there is M such that Xm and Ym contain both XL and Yq in their bounded interstice. But then by 



Lemma 4.5, the pair {Xm, Ym} is equal to {X' N , Y' N } for some N . To complete the proof that the 
algorithms eventually coincide, we need to show that we can choose M and N such that Xm = X' N 
and Y M = Y^. 

Suppose instead that Xm = Y' N and Ym = X' N . We may assume that Xm is smaller than Ym 
(otherwise, apply one more replacement to (Xm,Ym) to replace Xm with Ym+i and Ym with 
Xm+i)- It follows that Yjy is smaller than X' N , so that the next step of the algorithm, step (B), 
will set X' N+1 = Yjq and Y' N+l = X' N , and we therefore have Xm = X' N+l and Ym = Y' N+l . □ 



4.3. Self-similar packings. We also have already seen (Theorem 3.22) that V a is self-similar (but 
not the strip packing) if and only if a is quadratic over Q. It is a well-known fact that an irrational 
number a is quadratic over Q if and only if its continued fraction expansion is infinite and periodic; 
i.e., if and only if 

q = [ao, • • • , a„_i, cq, . . . , c m _i, cq, . . . , c m _i, . . .] = [ao, • • • , a n _i, Co, . . . , c m -\\. 



This fact, along with Theorem |4.7[ yields the following classification of self-similar half-plane pack- 
ings (which are not strip packings): 

Theorem 4.8. The similarity classes of self-similar half-plane (non-strip) circle packings corre- 
spond bijectively to finite non-repeating sequences of positive integers, up to cyclic permutation. 
(A sequence is non-repeating if it is not equal to the concatenation of multiple copies of the same 
smaller sequence.) 

Proof. Va is self-similar if and only if the continued fraction expansion for a is periodic, as we 
mentioned above. We identify the similarity class of V a with the minimal periodic part of this 
expansion. For example, if a = y/2 = [1,2,2,2,...] then we identify [V a ] with the one term 
sequence (2). If a = V3 = [1,1,2,1,2,1,2,...], we identify [V a ] with (1,2), or equivalently, 
(2, 1). Conversely, the finite non-repeating sequence (ao, . . • , a n -i) represents the quadratic number 
a = [ao, • • • , a n -i] , so that every such sequence represents a similarity class of self-similar half-plane 
packings (namely, the class including V a )- That each sequence represents exactly one similarity 



class follows from Theorem 4.7, since two periodic continued fractions expansions are eventually 



equal if and only if they have the same periodic part up to a cyclic permutation. □ 



Examples. Based on the classification in Theorem 4.8 we can give examples of the simplest self- 
similar half-plane packings. In some se nse, the simplest such packing is given by a = [1, 1, . . . ] = 

1+ 2 V ^ . From the perspective of Theorem 4.8 this is represented by the singleton sequence (1). The 



corresponding circle packing V a has a self-similarity constructed from a single circle replacement. 



This is displayed in Figure 11 



The next simplest example is a = [2,2, . . .] = 1 + v2, represented by the singleton sequence (2). 
The corresponding circle packing V a has a self-similarity obtained by doing two circle replacements. 



This appears in Figure 12 



There are two different self-similar packings for which a self-similarity involves three circle replace- 
ments. Corresponding to the sequence (3), we have a = [3, 3, . . .] = 3+ ^* - The packing V a is 



shown in Figure 13 Corresponding to the sequence (1,2) we have a = [1, 2, 1 2,...] = ^fi. This 



packing is shown in Figure [14} 

We have seen that the self-similarity of the packing V a is reflected in the continued fraction expan- 
sion of a. In fact, it turns out that every self-similarity of V a comes about from the periodicity of 
the continued fraction expansion. In particular, the self-similarity groups of the packings in Figures 
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^0 ^0 

Figure 11. Eight generations of the packing V a , where a = [1] = 1+ 2 V ^ . The circles 
corresponding to the circle replacement algorithm are shaded gray. 
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Figure 12. Eight generations of the packing V a , where a = [2] = 1 + y/2. The 
circles corresponding to the circle replacement algorithm are shaded gray. 
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Figure 13. Eight generations of the packing V a , where a = [3] = 3+ ^ 13 . The 
circles corresponding to the circle replacement algorithm are shaded gray. 
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TT]|T4 are generated by those arising from the circle replacement algorithm. We can make this more 
precise as follows. 

Suppose for simplicity that a has a purely periodic continued fraction expansion: 
a = [ao, • • • , a n _i, ao, • • • , On-li °0) • • • ] = [«o> ■ • • , a-n-i]- 



(By Theorem 4.7 we lose no generality in doing so.) After performing N = oq + • • • + a n -i + n 



steps of the continued fraction algorithm, the initial segment consisting of the periodic part of the 



continued fraction expansion is removed; therefore on = a. It follows by Lemma 4.2 that applying 
N steps of the circle replacement algorithm to the pair (Xq,Yq) = (Cn ), C(o,i)) yields a pair 
(Xn, Y/v) whose curvatures are in the same ratio (i.e., a 2 ) as the original circles (Xq, Yq). There is 



therefore a similarity of V a that maps Xq to Xjy and Yq to Yjy by Lemma 3.5 (Note that if the 
continued fraction expansion of a is not purely periodic, the above argument determines instead a 
similarity mapping (Xm, Ym) to (Xjy, Y/v) for some M, N.) In fact, this argument may be repeated 
to show that there is a similarity of V a that maps Xq to X^n and Yq to Y^n for each k G Z>o- The 
next lemma tells us that every self-similarity of V a arises from the continued fraction expansion of 
a in this way. 

Lemma 4.9. Suppose a is the positive real number with periodic continued fraction expansion 
[ao, • ■ • , a n _i] and that this is the minimal periodic part. Let <3? denote the generator o/Symm('P Q ,) = 
Z for which <I> has scale factor greater than 1, and let (Xq, Yq) = (C( lj0 ), C(o,i)) ^ e the generators 
for Vol ■ Then 

$ k (X ) = X kN , $ k (Y ) = Y kN 

for each k > 0, where N = oq + • • • + a n __i + n as above, and § k denotes the k-fold composition of 
$ with itself (and <£° is the identity map on M?). 



Proof. We have already argued that for each k > there is a self-similarity & k of V a such that 

<S> k (X ) = X kN and $ fc (Y ) = Y kN . 

Furthermore, if N does not divide m, then a m ^ a, so that the only self-similarities that correspond 
to the circle replacement algorithm are the <3?fc. 

Because the radius of Yq is 1, the scale factor \x k of § k is equal to the radius of ^ k {Yo) = Y k ^, which 
increases as k increases by the definition of the circle replacement algorithm. (Note that the Y k jy 
must necessarily be distinct.) Similarly, since the scale factor for <£ fc is fi k , with fj, > 1, the scale 
factor for <J? fc also increases as k increases. It will therefore suffice to show that each self-similarity 
<j> with scale factor greater than 1 satisfies 

4>{Xq) = X m and (p(Y ) = Y m 

for some m > 1. 



By Lemma 4.5 we may reduce this problem to showing that each of Xq and Yq either lies in 
the bounded interstice for 4>{Xq) and 4>(Yq) or is equal to one of 4>(Xq) and <j){Yo). Because self- 
similarities preserve the basic structure of the packing, Xq (resp. Yq) lies in the bounded interstice 
for <I>(Xq) and (^{Yq) if and only if 4>~ 1 (Xq) (resp. (/)~ 1 (Yq)) lies in the bounded interstice for Xq 
and Yo- It is therefore enough to show that if tp is a self-similarity of V a with scale factor less than 
1, then iJ;(Xq) lies in the bounded interstice for Xq and Yq. In that case i^{Yq) will have to either 
be in the bounded interstice, or be equal to Xq or Yq. 



By the proof of Theorem 3.22, for any nontrivial self-similarity ip of V a , we have 

ip(X Q ) = -0(C(i,o)) = C^-e^ - w ) 
ip(Y ) = V(C(o,i)) = C^w^ay 
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where x and y are integers satisfying \x 2 — Dy 2 \ = 4 with y 7^ 0. Here f(x) = px 2 + qx + r, p > 0, 
is the primitive integer polynomial satisfied by a, and D = q 2 — Apr is the discriminant of /. 

To show that ^(Cn i0 )) either lies in the bounded interstice for C(i )0 ) and C(o,i)) or is equal to C( ,i), 



it will suffice by Lemmas 2.15 and 2.20 to show that 



(4.10) x m > and _ yr> 0. 

To see this, we first recall a result due to Galois [HE] concerning purely periodic continued fractions. 
Because a has a purely periodic continued fraction expansion, it is a reduced quadratic number; 
i.e., a > 1 and — 1 < a' < 0, where a' is the quadratic conjugate of a. Since q = —p(a + a') and 
r = pad , and since p > 0, it follows that g, r < 0. It is now sufficient to show that x > and 
y > 0. 

First note that we cannot have both x, y < since then the label on ^(Cn )) would consist of two 



negative numbers, which is impossible by Proposition 2.21 We now show also that x and y cannot 
have different signs. 

Let fjf be the scale factor of ip which, by assumption, is less than 1. Since curv(C( ,i)) = 1 and 
curv(^(C (0il) )) = (ypa + 5±K£) 2 ; , it follows that 

(ypa + -?p J 

Therefore 



> 1. 



2 

Since a has purely periodic continued fraction expansion, it is greater than its conjugate, so we 
have a = ~ q ~2^f^ ■ The above inequality then implies 



\x + yVD\ > 2. 

Now x 2 — Dy 2 = ±4, so x = ±\JDy 2 db 4. If x and y have different signs, then this inequality 
becomes 

\^jDy 2 ±A-^fDy 2 \ > 2. 

But I \/t + 4 — -^1^2 for all i > 0, so in fact x and y must have the same sign (and y 7^ since 
we assumed that tp was a nontrivial self-similarity). This completes the proof. □ 

Finally, we can also use continued fractions to see which packings have orientation-reversing self- 
similarities. 

Theorem 4.11. For a positive real number a, the circle packing V a has an orientation-reversing 
self- similarity if and only if the continued fraction expansion of a has odd period. 



Proof. We know from Lemma 4.9 that all self-similarities of V a correspond to periods in the contin- 
ued fraction expansion of a. The self-similarity $ is orientation-reversing if the corresponding pairs 
of circles in the circle replacement algorithm, say (X n ,Y n ) and (X m ,Y m ), for which $(X n ) = X m 
and $(Y n ) = Y m , satisfy X n -< Y n , but Y m -< X m . Since the orientation of the circles X^^Y^ 
changes once for each case of step (B) in the algorithm, that is, for each term in the continued 
fraction expansion, we see that $ is orientation-reversing if and only if it corresponds to an odd 
period. □ 

In the examples, we see that the circle packings V a for a = [I], [2], [3] do have orientation-reversing 
self-similarities, whereas that for a = [1, 2] does not. 



APOLLONIAN CIRCLE PACKINGS OF THE HALF-PLANE 



35 



References 

1. H. Davenport, The Higher Arithmetic, Cambridge University Press, 2003. 

2. Evaristc Galois, Demonstration d'un theoreme sur les fractions continues periodiques, Ann. Math. Pures et Appl. 
19 (1829), 294-301. 

3. R. L. Graham, J. C. Lagarias, C. L. Mallows, A. R. Wilks, and C. H. Yan, Apollonian Circle Packings: Number 
Theory, J. Number Theory 100 (2003), 1-45. 

4. , Apollonian Circle Packings: Geometry and Group Theory I. The Apollonian Group, Discrete Comput. 

Geom. 34 (2005), 547-585. 

5. , Apollonian Circle Packings: Geometry and Group Theory II. Super-Apollonian Group and Integral Pack- 
ings, Discrete Comput. Geom. 35 (2006), 1-36. 

6. , Apollonian Circle Packings: Geometry and Group Theory III. Higher Dimensions, Discrete Comput. 

Geom. 35 (2006), 37-72. 

7. G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, sixth ed., Oxford University Press, 
2008. 

8. Michael Jacobson and Hugh Williams, Solving the Pell Equation, CMS Books in Mathematics, Springer, 2009. 

9. Daniel Pedoe, On a Theorem in Geometry, Amer. Math. Monthly 74 (1967), 627-640. 



